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Resolvents and Seiberg-Witten representation for 

Gaussian /3-ensemble 

A.Mironov* A.Morozovj A.Popolitovj Sh.Shakirov § 

ABSTRACT 

The exact free energy of matrix model always obeys the Seiberg-Witten (SW) equations on a complex curve defined by 
singularities of the quasiclassical resolvent. The role of SW differential is played by the exact one-point resolvent. We show that 
these properties are preserved in generalization of matrix models to beta-ensembles. However, since the integrability and Harer- 
£\J ' Zagier topological recursion are still unavailable for beta-ensembles, we need to rely upon the ordinary AMM/EO recursion to 

evaluate the first terms of the genus expansion. Consideration in this paper is restricted to the Gaussian model. 
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1 Introduction 

Seiberg-Witten (SW) prepotentials J- (a) [TJ [2j [3j H] are defined from the peculiar set of implicit equations: 



3 = fofi 



< (i) 

Here is an (m, 0) analytic form (holomorphic, meromorphic or even possessing essential singularities) on a 
J> , family of d = 2m complex manifolds with a system of conjugated (m, 0)-cycles A and B. When the system (JT]) 
is resolvable (its consistency is guaranteed by the Riemann identities) then a are called flat coordinates on the 
moduli space of the family (or simply the flat moduli), and J 7 (a) is a "quasiclassical" or Whitham r-function, 
on this space, satisfying a set of the (generalized) WDVV equations (usually as a consequence of the "residue 
formula") [5]. This is by now a classical branch of science, presented in big detail in numerous papers. 

A little more recently it has been realized that despite the "quasiclassical" nature of the SW equations, 
they perfectly survive various "quantization" procedures. The true conceptual meaning of this phenomenon 
still lacks understanding, but the very fact is getting established more and more reliably. The latest example 
is the Bohr-Sommerfeld representation [6l [7] of the LMNS free energy [8] in the Nekrasov-Shatshvili limit [9] 
£2 = 0: if also ei = 0, then this free energy is just the ordinary SW prepotential of [U|3], but remarkably eqs.([T]) 
survive when at least the first "quantization" parameter t\ is switched on. Actually it is claimed in [TOl E] that 
they will survive even further: when both t\ and ti are non- vanishing. And this claim is inspired by the AGT 
relations [T^J [T21 [TJ] , which provide a matrix model representation of the LMNS partition function [TJ1 [TU] . 
Then one can use the previous fact of the fundamental importance: that the exact matrix model free energies 
possess the SW representation with the role of SW differential played by the one-point resolvent 



Z 1X1 I MM 



which is a meromorphic differential on the spectral curve £ MM . Again, the SW representation is a kind of 
straightforward for the planar free energy (where it is actually discussed since |16[ 117] ) but, remarkably, it 
survives when all higher-genus corrections in powers of the string coupling constant gs (i.e. the t'Hooft's 
coupling A = gN) are switched on. This fact is still less known and under-appreciated. It was mentioned in 
passing in [TB] and in [TS] , but its real significance can be illustrated by the recent suggestion of [TT] to use it in 
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a conceptual proof of the AGT relations: by applying the topological recursion procedures |T8} to construct 
a double deformation of the original SW free energy, to g s = \/ — e i £ 2 ^ and to (3 = b 2 = —ei/e 2 ^ 1. 

It is the goal of the present paper to provide more illustrations to the SW representation of exact matrix 
model free energies and to make this crucially important technique more understandable and useable. The 
paper is dedicated entirely to this issue and we avoid mixing it with the other subjects. The first illustration of 
this kind was already provided in the Appendix to [TT], we reproduce it here and extend to non- unity j3. We do 
not address non-Gaussian models here, since this requires usage of rather heavy techniques, but we will address 
this question in forthcoming papers. Of course, for AGT applications one needs an essentially non-Gaussian 
/3-ensemble: the open-contour Dotsenko-Fateev integral a la [TU] [2TJ , which we do not consider in the present 
text. However, the SW representation undoubtedly exists there as well, also for arbitrary /3 and in all orders of 
the genus expansion. 

2 The case of (5 = 1: a source of questions and educated guesses 

The partition function is defined as 

Z{N) = -L [ dXx . . . dX N J] |Ai - X-j \ 2P e~h ^ A * (3) 

" J i<j 

For (3 = 1 it is equal to 

N-l 

Z(N) = ^ N V9 N ' 3 I] k - ( 4 ) 
fe=l 

Hence, for the free energy F = hiZ one has (ignoring the terms quadratic and linear in N) 

N-l 



F(N) = ]T ln(fc!) (5) 



fc=i 



It turns out [22] that Z(N) is a Toda-chain r-function and F(N) possesses the Seiberg-Witten representation 
That is, let pi(z) be the one-point resolvent of the model 

^) = (eA) ( fi ) 



Then the system of partial differential SW-equations 



~ TT~ / Pi{z)dz = a - (f p\{z)dz = , (7) 

2m J A J B da 

determines the SW prepotential, which, as one can check using the explicit expression for the resolvent from 
[TT] . is equal to the free energy 

F SW (N) = F(N) (8) 
and this equality just gives the SW-representation of free energy of the matrix model. 



Another remarkable fact is that the one-point resolvent satisfies the difference equation [26J [27] 

dz 

which implies that its B-periods satisfy [UJ 



Pl (N + l,z)+ Pl (N-l,z)- 2 Pl (N, z) = (N, z), (9) 



Tl B (N + l) + Tl B (N-l)-2Tl B (N) = ~^ (10) 

Eq.Q is closely related to integrability of Z(N), that is, to the Toda chain equation [22] 

Z(N)d 2 Z(N) - dxZ(N)diZ(N) = Z(N + 1)Z(N - 1), (11) 



where dxZ{N) = (£. A*) and dfZ(N) = ^ Xtf 
Eq. pT?]) is a weaker corollary of ©• 

Knowing these facts, the following questions arise naturally: 
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• Does © hold as well in the ft ^ 1 case? 

• Is there some /3-deformed version of © and (ITU|) ? 

The rest of the paper is devoted to the affirmative answer to the first question. A partial progress in 
answering the second one is outlined in the Appendix. 



3 Resolvents 

3.1 Ward identities: generalities 

A powerful technique for evaluating correlators in matrix models is known under different names: of the Virasoro 
constraints, of the loop equations, of the Ward identities [23, IS] . It relies on "the general covariance" of partition 
function: that is, the invariance of integral under arbitrary change of integration variables. For the eigenvalue 
model, not obligatory Gaussian, the Virasoro constraints can be deduced as follows [21] ■ Consider the obvious 
identity 

J2 fd\ 1 ...dX N -^-(xtA^e-^ v ^S il ...S im )=Q, (12) 

k 

where Si — J2 a Ki an d A is the absolute value of the Van- der- Monde determinant. Here V(X) — J2k^kX k ; in 
the Gaussian case only T2 = 1/2 is non- vanishing. 
One can easily check that 

J2 4- (AJJA 2 ' 9 ) = (p ]T S a 5 B _i_„ + (1 - P)nS n . x ] A 2/3 (13) 

k k V a=0 / 

and this is the only piece of equation that changes when one changes /3. 
Differentiation of the potential term gives 



(14) 



and this is the only model-dependent part of our consideration. 
Differentiation of the remaining terms gives 

d m 
^fc "TT~ (Sii ■ ■ ■ Sim ) — h'Sii ■ ■ ■ Sij+n-l ■ ■ ■ &i m (15) 
k ° Ak 3=1 

Now, having all the ingredients of the equations, one can write them in various forms. 

Virasoro constraints. If one denotes the disconnected correlator as 

C io ,., im ={S ia ...S im ) (16) 

the above considerations imply that 

n—l ^ rn 

ft 2^ Ca,n-l-a,ix,.,i m + (1 — 0)nCn—1,ii.,.,i m 5^ kTkC n -X+k,ix,.,i m + J]] ijCi^.^j+n-l,.,!,,, = (17) 

a=0 & k j=l 



Differential (W) operators. If one works with the generic partition function (with infinitely many non-fixed 
times) one can write these equations as a differential equation on the (full) partition function. Namely, let the 
potential have the form 

oo 

V(X) = (T a + t a )N + Y / (Tk+t k )X n , (18) 

fe=i 
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where are background values of source fields (usually, only finitely many of them are non-zero) and are 
perturbations of these background values. The partition function is thought of as a formal series in tk- Note 
that, for the non-normalized average, one has 

(S a ) = -g^-{l) (19) 
and, hence, the Virasoro constraints (fT7|) can be written as 

ei- ■ ■ ■ et- (J>< T > + "te + 9(1 - » + " £ sat:) z - (20) 

Loop equations. They equations arise when one sums up all the Virasoro constraints with the weights —^ft 
and writes the resulting equation in terms of the resolvents. For this purpose, it is convenient to rewrite the 
Van-der-Monde part of the identity as 



Ej^KA") - s^l^U^r' l a* (2i) 



Now summing up all the contributions one gets 

a i \ , /a -,n 5 / n , 5 r(zi,.,z m ) -r(zi,.,z ,-,^m) , 00 x 
Pr(z ,z ,zi, .,z rn ) + {p- l)—r(zo,zi, .,z m ) + V ^ (22) 

J=l J 

oo oo A; — 2 

kT k z^ 1 r(z , zi, ., z m ) + - ^ kT k ^o^ - 9 dzz k ~ 2 ~ 3 'r(z, Zi, ., z m ) = 0, 

where r(zo, ., z m ) is the disconnected resolvent 

To solve these equations perturbatively in g one has to rewrite the disconnected resolvents in terms of the 
connected ones. Then, the iteration procedure becomes well-defined: at each step of the procedure one has a 
system of linear equations for pi j, with fixed value of i + j. The expansion in powers k of g, as usual, counts 
contributions of genus k/2 Riemann surfaces in string (or topological) expansion. Here Pij stands for the genus 
j contribution to the i-point connected resolvent. 



3.2 Prerequisite: particular correlators 

The Ward identities in the form of the Virasoro constraints are very helpful in evaluating individual correlators 
Cii...i m - The advantage of this method is that the answers are exact in g and one may not rewrite the discon- 
nected correlators in terms of the connected ones for the iteration procedure to work (this simplifies the work 
drastically if one uses symbolic computer computations). 



To give an impression of what individual correlators look like we provide the first few one- and two-point 
correlators. Note that K denotes the connected correlators, and A = Ng. 



K k = ^ H£A fe ) = «£A fe )) (24) 

K Q = A 

K 2 = A(/3A - j8 + 1) 

K A = A (2/3 2 A 2 - 5/3 2 A + 3/3 2 + 5/3 A -5/3 + 3) 

K 6 = 5/3 3 A 4 + (22/3 2 - 22/3 3 ) A 3 + (32/3 3 - 54/3 2 + 32/3) A 2 + (-15/3 3 + 32/3 2 - 32/3 -I- 15) A 

K 8 = 14/3 4 A 5 + (93/3 3 - 93/3 4 ) A 4 + (2 34/3 4 - 398/3 3 + 234/3 2 ) A 3 

+ (-260^ 4 + 565/3 3 - 565/3 2 + 260/3) A 2 + (105/3 4 - 260/3 3 + 331/3 2 - 260^ + 105) A 
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Kkj = C kd - C k C 3 = ^Ei Ei m)) (25) 

#i,i = A 

#1,3 = 3A(/3(A- 1) + 1) 

#2,2 = 2A(/3(A- 1) + 1) 

Ki, 5 = 10/3 2 A 3 + 5 (5/3 - 5/3 2 ) A 2 + 5 (3/3 2 - 5,3 + 3) A 
#M = 4A(/3(A - 1)C#(2A - 3) + 5) + 3) 

tf 3 ,3 = 3A(/3(A - 1)(/3(4A - 5) + 9) + 5) 

In terms of the CFT-inspired variables M — bA and Q — b — ^, b = y/]3 these read 

= f (26) 

K 2 = M(M - Q) 

Ki = Mb(l + 2M 2 - 5MQ + 3Q 2 ) 

K 6 = Mb 2 (5M(2 + A/ 2 ) - (13 + 22M 2 )Q + 32MQ 2 - 15Q 3 ) 

# 8 = Mb 3 (21 + 14M 4 - 93A/ 3 Q + 160Q 2 + 105g 4 - 5AfQ(43 + 52Q 2 ) + Af 2 (70 + 234Q 2 )) 



= f (27) 

#1,3 = 3M (M - Q) 

# 2 , 2 = 2M(M - Q) 



K h5 = 5M6(1 + 2Af 2 - 5A/Q + 3Q 2 ) 
K 2A = 4Mb(l + 2M 2 -5MQ + 3Q 2 ) 
#3,3= 3Af6(l + 4Af 2 -9A/Q + 5Q 2 ) 

Note the remarkable simplification in comparison with (|24l) and (l2~5j) . 



3.3 The answer for resolvent at /3 — 1 

Just for completeness (and in part to emphasize the relative complexity of the (3^1 case) we begin from the 
well-known one-point resolvent at j3 = 1 |18j : 



The particular genus contributions are 



z — Xi , 

I k=0 



pi,o(z) = \ ( z -y( z )) 



Pi,2 (z) = ^7-5 (29) 



A 

y(z) 5 

2lA(A + z 2 ) 



Pia(z)= , 

2/(z) n 
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, , 11A(158A 2 + 558Az 2 + 135z 4 ) 
(*) = ^Tr • 

(30) 

where y(z) 2 — z 2 — 4A and all pi,2fc+i vanish. General formulae for pi ; 2n can be obtained by integral transfor- 
mation from exact Harer-Zagier functions, see [26l [T8l [27] . 

3.4 The answer for pi at generic (3^1 

The loop equations (|22|) in the case of Gaussian model acquire a very simple form 

P> (z , z , zi, 2m) + {p- 1)-^ — T (*0, Zl, ., z m ) + > J - (31) 

OZq OZj Zj — Zq 

1 A 

z r (Zq, Zi, ., Zm) + -~T (zi, ., z m ) = 

9 9 

where r denotes the disconnected resolvent 

r (z lt ., z m ) = (J2 ■ ■ ■ E ) ( 32 ) 

To solve this system of equations one should rewrite the disconnected correlators in terms of the connected 
ones and substitute the connected correlators by their Laurent expansion [25] . 

Thus, assuming that pi (z) — ^ X)i*lo Phi ( z ) ' 9 % ( so that the even parts of p are associated with oriented 
surfaces, while the odd parts with the non-oriented ones, with half-integer genera), one gets for the first few 
terms: 

z y(z) _ 1 



I L A £ a i 

2 2/9 2/3 2 _ p — 1 



j/(z) y(z) 2 2/3y(z) V y(z) 



. 5/3 2 A~9/3A + 5A g + | - 2 -fo - f + 2z 

Phi (2) = 1 7-^5 — H t^t (33) 

y(z) 5 ?/(z) 3 y(z) 4 



*..(*) = (/5-i)(^^^(i-^ T ) + 



5A(5-9/3+5ff 2 ) , -Az(30-43/3+30/3 2 ) 

Riy 1 Rp 



_ 37 /3 3_27!ei + i99 / 3 + 37 _273 _ 37 ,3 3 z + HM^ _ 199(32; _ 37* + 273* 419/3 4 A- 1357/3 3 A+ 1897/3 2 A- 1357/3A+419A 
W - RP + RP + R^F 

-240ff 4 Az+824ff 3 Az-1168/3 2 Az+824/3Az-240Az , 1105fi 5 A 2 -3240/3 4 A 2 +4375/3 3 A 2 -3240,3 2 A 2 + 1105/3A 2 
"•" Rip + y(z)H 



Pi, 5 (z) = (/?-!) 



(706-2379/3+3367/3 2 -2379/3 3 + 706/3 4 ) , 4351-13458/3 + 18508ff 2 -13458/3 3 +4351ff 4 
2/3y9(z)(l-^ y ) + 



3Az(1530-4241/3+5764,3 2 -4241/3 3 + 1530,3 4 ) , 55/3A 2 (221-648,3+875fi 2 -648ff 3 +221/3 4 ) 
Rip 1 Rip 

4/3 A 2 z(3390-7883/3+ 10420/3 2 - 78830 3 +3390,fl 4 ) " 
RIP 

4081/3 5 - 404 ° 5 ' 34 H-44699/3 3 -57155/3 2 +44699/3+ ^f- - ^jg -4081/3 5 z+ 4040 2 5g4 ' -44699/3 3 z+57155/3 2 z-44699/3z-^|^ + fflf^ 
Pi, 6 (2J — ^fjp I " 7/(Jp 

i 77597/3 6 A-340402/3 5 A+702694/3 4 A-878293/3 3 A+702694/3 2 A~- 340402^ A+77597A 

Rip 

i -59040/3 6 Az+269328/3 5 Az-564000/3 4 Az+707424£i 3 Az-564000/3 2 Az+269328/3Az-59040Az 

+ RIP 

, 451720/3 7 A 2 -1792898^ 6 A 2 +3483419/3 5 A 2 -4266464^ 4 A 2 +3483419/3 3 A 2 -1792898/3 2 A 2 +451720^A 2 
+ RIP 



~189840/3 7 A 2 z+821128/3 6 A 2 z-1656256/3 5 A 2 z+2049936/3 4 A 2 z-1656256,3 3 A 2 z+821128,3 2 A 2 z-189840/jA 2 z 



828250/3 8 A 3 -3012930/3 7 A 3 +5531740/3 6 A 3 -6644070/3 5 A 3 +5531740/3 4 A 3 -3012930/3 3 A 3 +828250ff 2 A 3 

RIP : 
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here y (z) — z 2 — 4A/3 defines the spectral curve, which in this case is the torus with a degenerated handle 
(located at infinity of the complex plane) . 

(3 —> 4 symmetry. The AGT relation implies that the /3-deformed matrix model should be related to some 
CFT, with the central charge of the corresponding CFT given by 

c=l- 6 (^-i) 2 (34) 

This hints that there should be the symmetry (3 — > -g present in the matrix model despite this is far from obvious 
in the original expression (3). And, indeed, one can see that if one rescales the quantities in the following way: 

z' = V/?z, p'i, g = \fp 9+ p\, g , (35) 
the resulting p[ g are symmetric w.r.t. (3 — > -e. 

4 Seiberg-Witten construction 

4.1 Ideology 

The Seiberg-Witten construction, originally proposed to obtain the low-energy effective action in M — 2 SUSY 
gauge theory is in fact a manifestation of a more general statement. 

The starting objects in the SW representation arc the algebraic curve and the meromorphic differential Xsw 
on it. Given such a data, one writes the following system of equations 

X sw ~ flt / Asw ~ i (36) 

where Ai and Bi form a symplectic basis of cycles on the algebraic curve and proportionality coefficients in 
equations slightly depends on setting. 

It turns out that a huge source of the SW data is provided by the eigenvalue models (EVM). Namely, the 
algebraic curve is the spectral curve of the given EVM, while the SW differential is pi(z)dz, p\ being the one- 
point resolvent. Note that the original SW construction corresponds to the zeroth order of genus expansion 
of the resolvent in <?, and taking into account further terms of the expansion corresponds to the deformation 
(quantization) of the original SW differential and prepotential. Remarkably, the all genus free energy, not only 
its genus zero part, continues to satisfy the SW equations. 

We fix the proportionality coefficients in the SW equations as follows 

- 77— / Pi(z)dz = ^ - (3 (£ pi{z)dz = d ^ sw | (37) 

27T? J A . J B . dai 

as relation between the free energy and the SW prepotential is most transparent in this way. Note that j3 
appeared as a coefficient in the second equation, [7]. 

4.2 Calculation of A- and B-periods 

Now let us apply the SW construction to p\ that we found in section I3~4l 
The spectral curve is given by the equation 

y 2 = z 2 - 4A^ (38) 

The A-cycle encircles the ramification points — ^/4A/3 and y/4A/3, while the B-cycle encircles y/4A/3 and 00. 

Since in this case the value of A-period is equal to the residue at infinity, the A-pcriod gets contributions 
only from pi t Q and pi t \: 

1 1 - 8 

a = -—f p{z)dz = N+—f (39) 

Note at this point that the dependence a(N) is linear and one can safely substitute -g- by -J^ in the SW 
equation to simplify calculations. 
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Evaluating the B-periods is more tricky, the following formula is of great use 



+ °° dz i r(p - i)r(i - p/2) 



(40) 



To deduce this formula, one has to make the change of variables z — V4A/3 and notice that the resulting 



integral is proportional to the integral representation for the Euler i?-function 

n + OO /> + OC 



(z 2 - 4/3 A) - p/2 dz = (4/3Ay p/2+1/2 J> (w 2 - l) 



1-P/2 



dw 



(41) 



= (4/3A)" p/2+1/2 4-P/ 2 2 /\l - Cr p/2 C p ' 2 dC = 2- 2p+3 (/3A)- p/2+1/2 F (l 2 ] T p . iP 1] 

The terms in (|33|) with odd powers of z do not contribute to the periods, since they are total derivatives. For 
instance, 



zdz 



1 



d 



1 



= 



Note that this is a contour integral and the contour does not pass through the singularities of the integrand 



(42) 



Alternatively, one may exploit the fact that y(z) satisfies the differential equation 

■^f = -2/W" 2 , 

and so do its B-periods. Together with the initial conditions 

dz 

ib y{z) 

and 



= -ln^A 



y p dz = 



A=0 



this gives (only minor modifications occur in comparison with (3 = 1 case) 
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-2/3(A-Alog(/3A)) 



-10/3 



6/3 /3A 2 log(/3A) 



11/TA J 



2^ 2 A 3 log(/3A) 



-14/3 (5/3 3 A 4 log(/3A) - ^g^-j 



18/3 f 959 ^/" - 14/3 4 A 5 log(/3A) 



22/3 42/3 5 A 6 log(/3A) - 



.26/3 ( 1197 3 9 f A? - 132/3 6 A 7 log(/3A) 



1^ 



-log(/3A) 



i 

'2/3A 



12/PA 5 



60,3 3 A 3 



280/3 4 A 4 



1260,3 5 A 5 



5544£ 6 A 6 



(43) 

(44) 
(45) 



Thus, for the B-periods of p\j one gets (in the case of and p\ t \ one has to evaluate the integrals for 
p = — 1 + e and 1 + e respectively, and then to neglect the terms which diverge as e — > 0; since these terms are 
constant and linear in A, this is safe) 

Pi,o{z)dz = —A In A 



(46) 



pi,i(z)dz = i— -In A 



pi !2 (z)dz 



pi, 3 (z)dz 



-1 + 3/3-ff 2 
12/3 2 A 

1-/? 
24/3 2 A 2 
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1 + ^4 

b ' 240/3 4 A 4 

2 + 10 1 + 7f3 4 - 2f3 6 



pifi(z)dz = 



2520/3 6 A 5 



Remarkably, despite the complexity of p% i grows very fast (exponentially) with increasing i, their B-periods 
complexity increases not so fast (linearly). 

Already at this stage one can see that these formulas agree with the generic ones from [25] 



m+l 

Pl,2m+2 = 



V> T(2m + 1) 2s 1 



s=0 



- i J_ J—\ B 2m+2 (2m-l) 1 
Pi,am+i-l 2/3 2/?2m j ( 2TO+ l)(2m + 2)7V 2 ™' " ' [ ' > 



IB 

In |25j they were deduced from eq. (|5ip , which we are now aiming to derive. 
4.3 Relation to free energy 

Partition function for the /3-deformed Gaussian eigenvalue model is defined as 



z{n) = -LJd\ 1 ...d\ N Y[ l A * - a j \ 2P e ~* El x ' ( 48 ) 



and can be calculated explicitly. Generalization of (4) for fj ^ 1 is (see 

AT 

-N 



z(n) = V2V V jgW+a-v" rr £0 + W . i_ (49) 

v y vy ii r(i + /3) r(iv + i) v 7 



Now we are ready to check that the free energy 

N 



F(N) =lnZ~^lnr(l + /3/c)-lnAn, (50) 



fe=i 



is equal to the SW prepotential. 

Indeed, one can calculate the A^-derivative of F(N) and apply the Euler-Maclaurin formula, see eq. ([72"T) in 
the Appendix, to obtain 



\ 7/ * ~ 12/3A 360FA 3 



g_„/A 



„4 l-j3 3 , „5 2-7/3 2 -7g 4 +2/3 6 , / IN 
» 240£ a A 4 ^ » 252O0 5 A 5 r u \iPi 



this expression can be now compared with (45), taking into account the factor —/3 in ((37] 
Finally, one obtains 



F = F sw (51) 



This is the main statement of the paper: 

The e xact free energy of the Gaussian /3-ensemble satisfies the SW equa- 
tions (1371) with the exact resolvents in the role of the SW differential. 
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A Appendix. Towards understanding of /3 ^ 1 

In this Appendix we outline a few topics which are poorly understood but are of crucial importance for the 
future theory of /3-ensembles. 



A.l Integrability 

At (3 = 1 we saw that the free energy and the resolvent satisfy integrable differential-difference equations © 
and (1101) . These equations are intimately related with the Toda integrable structure of the Gaussian matrix 
model (in general case, it becomes KP integrability). In particular, the Toda equation can be written as follows 



dt\ y ' Z 2 {N) 



2 In Z(N) = -i >- (52) 



and, in terms of the free energy, 



F(N + 1) - 2F(N) + F(N - 1) = In ( — ^ F(N) ) (53) 



Then, by differentiating w.r.t. t{ and applying the Virasoro constraints, one obtains 

Ki(N + 1) - 2K l (N) + Ki{N - 1) = t -^^±K i _ 2 {N) (54) 

Summing these equations with the weights —ptt, one obtains eq.((9j). Eq. ([T0| then follows from ([9]) if one 
integrates it along the B-period on the spectral curve. 



What is the /3-deformation of Toda/KP integrability is a very important and intriguing question, but hard 
to tackle straightforwardly. As we shall see, eq.([T0|) can be /3-deformed but integrability requires more: a 
/3-deformation of Q which is still unknown. 



A. 1.1 Difference equation for periods 

For (3 = 1 the equation (flU)) reads 

n B (A + i)-2n B (A) + n B (A-i) = -i (55) 

where II B (A) stands for the B period of p. 

Experimentally one can see (e.g. expanding the l.h.s. into the 4 series) that, for /3 7^ 1, this equation 
deforms to 

n B ( A + ^) - n * ( A ) - n B (a + + n B (A - 1) = --L ( 5 e) 



A. 1.2 Difference equation for resolvents 

We, however, were unable to find a /3-deformed version of ([9]), even the corrections of the first order in /3 — 1 are 
missing. What one can say for sure is that in the required generalization the both sides of ^ deform stronger 
than they do in (|56[) . 



10 



A. 2 Harer-Zagier topological recursion 

A detailed description of the Harer-Zagier functions for ft = 1 can be found e.g. in [57] ■ Description of matrix 
model correlators in terms of the resolvents has two advantages: it provides the Ward identities (|T7|) in a simple 
form of the loop equations (|22p and it reveals the important hidden structure, the spectral curve. The drawbacks 
are the divergency of series for the genus expansion and the lack of explicit formulas for exact correlators. 

The last two problems are resolved, e.g., by switching from the exact resolvents to the Harer-Zagier functions, 
where the correlators are weighted with additional factorial factors, i.e. by summing up the series expansion by 
the Pade method. 

For ft 1 much less is known. So far we were able to obtain the Harer-Zagier functions only for specific 
values of ft ^ 1. Attempts to evaluate, at least, the hrst ft — 1 correction lead to some generalizations of the 
hypergeometric equations which is a hint that something conceptual needs to be done for the results to become 
simple for arbitrary ft. Below our preliminary results are summarized. 

The one-point Harer-Zagier generating function is defined as 



4/3 



oo oo 



k=0N=0 



N 



T - 1 



(2k - 1)!! \t + 1 



(57) 



where Cu [^,ft^ is the one-point correlator in the /3-deformed matrix model with matrix size equal to N/ft. 
The case of ft = 1. 



<t>(P 



1 



(58) 



This is the classical result by J.Harer and D.Zagier. 

It satisfies the differential equation derived from the integrability conditions 



(59) 



where tp 



and A 



r-l 
r+1 ■ 



Now it turns out that the two- and three-point Harer-Zagier functions can be found as well, and they are 
expressed through the arctangent function [27], i.e. remain elementary functions. 



The case of ft = 2. The Harer-Zagier function for SO(N) matrix model has the form 

t yi(r + i) arctan (r^ Ii - 



z — 1) \/t — z — z 



2,z,t) = 



1+(2-3t)z+(2-2t+2t 2 )z s 



r — z — Z*t 



and it satisfies 



2z + 2t 2 z + t 



2z 



+ {*■ 



T 2 Z 



d 



<t>(J3 



dz 

1 + (t-1)z + 



2,z,t) 



(t — z — z 2 t) 3 / 2 



t 2 + 2t + 2t 2 z + t 2 
~2z + 2(2z + 1)(1 - zt) 



with the initial conditions <p(ft = 2, z, r) 
The case of ft = 1/2. The Harer-Zagier function for the Sp(N) matrix model has the form 

( : 



cl>(ft = 1/2, z,t) = 



arctan 



2\/ z-\-rz\/l — zr 



2-z-2tz 



1 + T 



and it satisfies 







[(:■ 









oz 



(1 - Tzfl 2 

(1 + t)(2-tz) d 



dr 



c(>(ft = 1/2, z,t)=- 



(60) 



(61) 



(62) 



(63) 



with the initial conditions <j)(ft = 1/2, z, r) = 1 + (r + l/2)z + . . . 

One can see that in both cases of ft = 2, 1/2, which correspond to classical groups, the Harer-Zagier functions 
remain expressed in terms of arctangents. This, however, is not the case in the general situation. 
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The case of j3 = 3. The Harer-Zagier function for /? — 3 satisfies the differential equation 



(1 + 8z 2 t + 24z 2 r 2 + 9z 3 t - zt - 6z - 33z 2 )cf> + (llz 2 r - 18z 3 - 2z + 9z 4 t)|§ (64) 
+ (9z - 12z 2 r + 12z 2 t 3 - 9zt 2 )§^ = 1 - 4z - 4zt 

at particular values of z it becomes the hypergeometric equation and, hence, has no solutions expressed in 
elementary functions. So, presumably, what one is searching for is some clever deformation of the arctangent 
function from the previously described cases. 

We observe that as we move further and further away from /3 = 1, the complexity of results increases. 
Further work is needed to clarify the situation. 

A. 3 Identities for free energy 

It turns out that for /3 ^ 1 the Gaussian free energy has more structure than one could expect. 
A. 3.1 Definitions 

Let us define the partition function without -Ay factor. To avoid an ambiguity, let us denote all the quantities 
in this normalization with tildes. 

The partition function for the Gaussian model we are considering is 

Z(N, P)= J d\ x . . . d\ N Y[ (\i - Xj) 20 ^* A ^ = N\Z(N, 0) (65) 

i<j 

Instead of ([49|) we now have 

N 

2 r ( i + 

Defining the free energy as 

N 

F(N,0) =\nZ~ ^lnr(l + /3jfe), (67) 
k—l 

where the equivalence means equality up to terms quadratic and linear in the matrix size N (they can be 
absorbed into redefinition of (3 and g). 

A. 3. 2 Difference equation 

Thus defined free energy satisfies a certain difference equation. Consider 

G{N,P) =F(N,j3)-F(N -l,/3) = \nT(l + pN) (68) 



w n = v^ N v-r 2+{i - nN n T -^§ (66) 



then it is obvious that 



which implies that 



G(N,0)-G[N-^,f] =hOSJV) (69) 
^S (N , m -JLG(N-l,e) = L ,70) 



A. 3. 3 Exact relation between Tsw and F 

Comparison of (f5"6"j) and ((70)) gives 



F\N-^,0j = lnr(l + /3fc), 



(71) 



so the only peculiarity is in change of upper limit of summation. In the case of j3 = 1, it becomes N — 1 
and acquires a clear physical meaning: division of partition function by N\ implies that the eigenvalues are 
indistinguishable bosons. 
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A. 3. 4 Direct comparison of series 

However, it is still instructive not to appeal to this difference equation argument, but to act straightforwardly 
and look directly at the perturbative expansions at large N in order to see, what one can deduce from them. 
One way to obtain these expansions is to use the Euler-Maclaurin formula. 

Euler-Maclaurin formula. We need this formula in the following form: 

— J2 /(*) = /(AO - ^f'(N) + -f"(N) WQ f""{N) -...= £ -J d™f(N), (72) 

k m— 

where B m are the Bernoulli numbers, %frt m ~ srzr- ^he ^ ow limit in the sum is inessential, as long as it 
does not depend on N . In the following examples it is chosen to be k = 0. 
The first examples are 

/(*) = i m N = 1 

/(*) = * &^ = *-i 

f{k)=k 2 ^ NW-WN-!) = N 2_ N+ 1 ^ N 2_2_N + 2_ 

f(k) = k* ^ N2(N 4 - 1)2 =N*-lN 2 + ±N =iV 3 -M^ + ^+0 

Different series as they are. Here all equalities are considered up to terms linear and constant in N or A. 
Summing up contributions from different genera, one gets 

„4 1-P 3 , n 5 2~7/3 2 -7l3 4 +2l3 6 , I 1 \ 
240£ a A 4 y 2520,3 5 A 5 

Expanding F at various points one gets 

F f A 1 \ - MA In A I P- 1 In A I „ 1 -3/ 3 +/ 32 i „2 ,3-1 , 3 -l+5j3 2 -g 4 , 

+■9 240FA 4 " + 9 2520/3 s A B h ( A^J 

pYM - 1 ft A In A I 1+13 In A I g i+W+P 2 n 2 1+13 . „3 l-5/3 2 +/3 4 /-^s 

VS / -/?AlnA+ — lnA + g 12fjA g - g 360/3 i A :i + U 4 ) 

„4 l+<3 3 _ , 5 2-7/3 2 -7/3 4 +2,8 6 , M\ 
+■9 240/3 3 A 4 + 9 2520/3 s A B h ( A^J 

F /^ A 1^ - 1 RK In A 4- l ~P In A 4- „I=M±£! + n 2 -±^L 3 l-5/3 2 +g 4 , 

„4 g 3 -l , 5 2-7/3 2 -7/3 4 +2/3 6 , M\ 
+■9 240FA 4 " + 3 2520/3 B A B h ( A^J 

Interpretation. From the above series, one can easily guess the following relations 



d_~fA 1 
dN J 

_d_ / A 

dN J 

_d_ ^ /A 

9Ar J 



^)-*(H>)--*(-7-^)--*(7--' 



(75) 



The first equality is expected to hold from our previous difference equation analysis. 

It turns out that F with shifted arguments also satisfies the same difference equation. Indeed, 



- F (N,-0) + F {N - 1,-/3) = -lnT(l - /3N) , 

-lnr(-/3A r )+l n r(l-^A r ) = ln(- /3N) (76) 

For F(— N — 1, /?) one has to assume that the lower limit of summation is less than — N — 1 (which is rather 
weird) 

- F (-N - 1, P) + F (-N, P) = - lnT(l - 0N) (77) 

Normally, shifting the expansion point for some function does not lead to series similar to the initial one, but 
produces something which looks completely different. The fact that this is not the case may be an indication 
that some not yet discovered mathematical structure is present here. Perhaps, it is a peculiar feature of the 
Gaussian potential or, may be such equalities have more general character. It is interesting to see which of these 
unexpected identities survive generalization to non-Gaussian eigenvalue models. 
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